function cp = specheat(Species,T)

R = 8.314;

O2h = [3.28253784E+00;1.48308754E-03;-7.57966669E-07;2.09470555E-10;-2.16717794E-14;-1.08845772E+03;5.45323129E+00];
O2l = [3.78245636E+00;-2.99673416E-03;9.84730201E-06;-9.68129509E-09;3.24372837E-12;-1.06394356E+03;3.65767573E+00];

H2Oh = [3.03399249E+00;2.17691804E-03;-1.64072518E-07;-9.70419870E-11;1.68200992E-14;-3.00042971E+04;4.96677010E+00];
H2Ol = [4.19864056E+00;-2.03643410E-03;6.52040211E-06;-5.48797062E-09;1.77197817E-12;-3.02937267E+04;-8.49032208E-01];

CO2h = [3.85746029E+00;4.41437026E-03;-2.21481404E-06;5.23490188E-10;-4.72084164E-14;-4.87591660E+04;2.27163806E+00];
CO2l = [2.35677352E+00;8.98459677E-03;-7.12356269E-06;2.45919022E-09;-1.43699548E-13;-4.83719697E+04;9.90105222E+00];

N2h = [0.02926640E+02;0.14879768E-02;-0.05684760E-05;0.10097038E-09;-0.06753351E-13;-0.09227977E+04;0.05980528E+02];
N2l = [0.03298677E+02;0.14082404E-02;-0.03963222E-04;0.05641515E-07;-0.02444854E-10;-0.10208999E+04;0.03950372E+02];

C8H16O2h = [1.55349351E+01;5.99912854E-02;-2.58415306E-05;4.78465321E-09;-3.27482927E-13;-7.60581758E+04;-4.61023856E+01];
C8H16O2l = [-6.91263883E-01;1.03659150E-01;-6.99110558E-05;2.45512669E-08;-3.65222106E-12;-7.12346831E+04;3.86102343E+01];

C8H10h = [1.09952650E+01;3.94471414E-02;-1.60141065E-05;3.02159851E-09;-2.15018598E-13;-4.33350070E+03;-3.36830297E+01];
C8H10l = [1.41688302E+00;4.24899260E-02;2.90259016E-05;-5.79774769E-08;2.22785462E-11;-5.38782056E+02;2.16289435E+01];

if T > 1000
    crit = 'h';
else
    crit = 'l';
end

Data_name = [Species,crit];

Data = eval(Data_name);

cp = R*(Data(1) + Data(2)*T + Data(3)*T^2 + Data(4)*T^3 + Data(5)*T^4);

end

